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The analytic result for the S-matrix in the saturation regime including the running coupling is 
obtained. To get this result we solve the Balitsky and Kovchegov-Weigert evolution equations in 
the saturation regime, which include running coupling corrections. We study also the effect of rare 
fluctuations on top of the running coupling. We find that the rare fluctuations are less important 
in the running coupling case as compared to the fixed coupling case. 
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O ' I- INTRODUCTION 

o ■ 

. The Balitsky-Kovchegov (BK) equation [J, Q is a non-linear evolution equation which describes the high energy 
Qh' scattering of a qg dipole on a target in the case of fixed coupling. An analytic solution to the BK equation in the 
, saturation region has been found by Levin and Tuchin [sl . The BK equation can be viewed as a mean field version of 
more complete equation [l| where the higher correlations are neglected: The S'-matrix of the scattering of two QCD 
dipoles on a target is replaced in the BK equation by the product of the S'-matrices of the individual dipoles. Such a 
^""^ replacement is legitimate only in the absence of fluctuations in the light cone wavefunction of the target |4|. However, 
' in Ref . [5] was shown that rare fluctuations do change the result for the S'-matrix in the saturation region. 

Recently, the evolution equations which include running coupling effects have been derived by Balitsky and 
O i' Kovchegov-Weigert 0, They found that the running coupling corrections are included in the BK kernel by 
JL , replacing the fixed coupling as in it with a "triumvirate" of the running couplings. A more complete evolution 
^ equation has been studied by Albacete and Kovchegov [l3|, they have calculated in addition to the Balitsky and 
(~| ] Kovchegov-Weigert equations also the so called subtraction contributions. A numerical solution of the more complete 
I evolution equations were given in [l3l | . 
^ In this work, we will analytically solve these equations in the saturation region and obtain an analytic result for the 

<^ . S'-matrix. We find that the running coupling corrections modify the S'-matrix a lot as compared to the fixed coupling 
\^ ' case. Moreover, we study the effect of the rare fluctuations on top of the running coupling in the way as it was done 
\^ in Ref. fs'] for the fixed coupling case. We find that the rare fiuctuations are less important in the running coupling 
■ case as compared to the fixed coupling case. 

(N ; 

a\ ■ 

O ■ II- FIXED COUPLING CASE 

00 : 

; The BK equation P, Q gives the evolution with rapidity Y — ln(l/x) of the scattering amplitude S{x±^,yi_,Y) of 
^ ■ a gg dipole with a target which may be another dipole, a hadron or a nucleus. The BK equation is a simple equation 
to deal with the onset of unitarity and to study parton saturation phenomena at high energies. The analytic solution 
to the fixed coupling BK equation for the S-matrix deep in the saturation regime has been derived by Levin and 
c3 Tuchinjl]. This solution agrees with the one derived by solving the BK equation in the small S limit Q. In this 
section we will give a simple derivation of the BK equation and its solution in the saturation regime. 
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A. The BK equation 



In the high-energy scattering of a quark- antiquark dipole on a target, it is convenient to view the scattering process 
in a frame where the dipole is moving along the negative z-axis and the target is moving along the positive z-axis. 
Further we assume that almost all of the rapidity of the scattering, Y, is taken by the target. We denote the scattering 
amplitude of a dipole, consisting of a quark at transverse coordinate x± and an antiquark at transverse coordinate 
y±, scattering on a target by S{x±,y±,Y). Now suppose we increase Y by a small amount dY. We wish to know 
how S{x±,y±,Y) changes with the small amount dY. If the rapidity of the dipole is increased while that of the 
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target is kept fixed, then the dipole has a probability to emit a gluon due to the change dY. We now calculate the 
probability for producing this quark-antiquark-gluon state. In the large limit the quark-antiquark-gluon state can 
be viewed as a system of two dipoles - one of the dipoles consists of the initial quark and the antiquark part of the 
gluon while the other dipole is given by the quark part of the gluon and the initial antiquark. Using the dipole model 
the probability for producing the quark-antiquark-gluon state from the initial quark-antiquark state is 



gipoie 



dP = —^d z^dY- — , (1) 

where z± is the transverse coordinate of the emitted gluon. The change in the S'-matrix, dS, for a dipole-hadron 
scattering is given by multiplying the probability dP with the S'-matrix 

■SpS{xj_-y±,Y) = [ d^z^- ''^'^x,/^'* ^ S^^Hx±- zj_,z±-y±,Y)- S{x±-yj_,Y) , (2) 

where S'''^\x±_ — z±, zj_ — y±,Y) stands for a simultaneous scattering of the two produced dipoles on the target (see 
the first diagram on r.h.s of Fig. [T]). The last term in ^ describes the scattering of a single dipole on the target 
because the gluon is not in the wavefunction of the dipole at the time of the scattering (see the last two diagrams in 
Fig.[T]). 

It is hard to directly use Eq. ([2]) to study problems of parton evolution and parton saturation phenomena at high 
density and high energy QCD, since S^^^ is not known. Using the mean field approximation for the gluonic fields in 
the target 

S^^Hx±-z±,z^--y^,Y)^Six^-z^,Y)Siz^-y^,Y) , (3) 
one gets the Kovchegov equation 

^S{x^ ^y^,Y)^^ [ d^z^- ''-^^-jf^ ^ [S{x^ - zx, Y)S{z^ - y^,Y) - S{x^ - y^,Y)] . (4) 

dY 2tt^ J [x^ - z±_y{z^~ y±_y 

With N{x±^ — y±, Y) = 1 — S{x± — y±, Y), another useful version of the Kovchegov equation is obtained 

^N{x^~y^,Y)^^ [ d'z^- -^[N{x^-z^,Y) + N{z^~y^,Y) 

dY 27r^ J {xj_ - zj_)''(z^ - y^y 

-N{x^ - y^,Y) - N{x^ - z^,Y)N{z^ - y^,Y)] . (5) 

Eq. ([5]) has the following probabilistic interpretation: when evolved in rapidity, the initial quark-antiquark dipole of size 
x±—y± decays into two dipoles of size x± — zj_ and zj_ — y± with the decay probability {x±—y±y/{{x± — z±y{z±—y±y) 
which is usually called as BFKL kernel. These two dipoles then interact with the target. The non-linear term takes 
into account a simultaneous interaction of two produced dipoles with the target. The non- linear term prevents the 
amplitude from growing boundlessly with rapidity and ensures the unitarity of the scattering amplitude. When the 
scattering is weak iV — > 0, the nonlinear term N{xj_ — z \ . Y)N(z \ — y±,Y) can be dropped and the linear equation 
remaining is the dipole version 0] of the BFKL equation [III . 




FIG. 1: Diagrams corresponding to terms in the evolution equation 
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B. Solution to the BK equation in the saturation regime 

In the high-energy regime where unitarity corrections become important or S{x± —y±,Y) is smah, Eq. (j4]) is easier 
to use since the quadratic term S{x± — zj_, Y)S{z± — y±,Y) can be neglected in which case one needs only keep the 
second term on the r.h.s of Q giving 



_d_ 
W 



d z± 



{x± - - y±y 



:S{x± - y±,Y) 



(6) 



In the above equation, we have assumed that S is small which holds only when the dipole size is large compared to 
1/Qs- Therefore the lower bound of integration in ^ should restrict to the regime {x± — y±f' ^ l/Qs '^^^ to 
the regime {x± — z\_)'^ ^ (-^-L ^ ^ 1/Qs- the logarithmic regime of integration one gets 



— S{x^-y^,Y) = -2^n 



{x±-V±) 



1/Ql 



d{z±_ - y±_f 



rS{xi_ -yA_,Y) 



(7) 



Note that the factor 2 in the above equation comes from the symmetry of the two regions dominating the integral, 
either from l/Qs ^ f 1 <C r2 ~ r or l/Qs ■^r2 -^r, ri ~ r. Now it is easy to get the solution to Eq. ([7]) 



S{x±_ - yA_,Y) = exp 



c / aNf 
'2 



TT 



(Y-Yo) 



Sixi_ - y±,Yo), 



(8) 



where we have used 



g|(y)==exp 



(9) 



and 



QUYo)ix±~y±)' 



1. 



(10) 



Eq. ^ gives the standard result given in the literature [8j]. We have gone through such a detailed "derivation" of 
([8]) since one of the main purposes of the present paper is to show how Eq. ([8|) is modified once running coupling 
effects and rare fluctuations are included. 



III. RUNNING COUPLING CASE 



The BK equation only considers the resummation of leading logarithmic (LL) as XBj) corrections with a fixed 
coupling constant The running coupling corrections due to fermion (quark) bubble diagrams, which would bring 
in a factor of a^A^/, modify the evolution equation, which is not leading logarithms anymore. Once including UgNf 
corrections, the obtained contributions have to divided into two parts, running coupling part and the "subtraction" 
part. The first part has a form as the leading order BK kernel but with the running coupling replacing the fixed 
coupling and the second part brings in new structures into the evolution equation. 



A. Balitsky and Kovchegov-Weigert equations 

The evolution equation including higher order corrections reads [isj 

dS{x^-y^,Y) 



dY 



n[s] - s[s] . (11) 



The first term in r.h.s of (jlip . TZ, which is referred to as the 'running coupling' contribution and resums all power 
of oigNf corrections to the evolution. The TZ. has a form as the leading order one but with modified kernel which 
includes all effects of the running coupling 

n[S{x^-y^,Y)]^ J d^z^Kix±,y^,z^)[S{x^-z±Y)S{z^^y^,Y)-S{x^^y±,Y)] . (12) 
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FIG. 2: The higher order diagrams contribution to BK evolution. 



The BK kernel is modified because the propagator of the emitted gluon in the original parent dipole is now dressed 
with quark loops in contrast to leading order or fixed coupling one. This modifies the emission probability of the 
gluon but doesn't change the leading order interaction terms (see Fig.2A). 
Using N{x± — y±,Y) = 1 — S{x± — y±,Y), another useful version of (fT2l) is: 



n[N{x^~y^,Y)] 



(fz^ K{x^,y^, z^) [N{x^ - z^,Y)+N{z^- y^,Y) 
- Nix^ - y^,Y) - N{x^ - z^,Y)N{z^ - y^^Y)] 



(13) 



with modified kernel K{x±^y±, z±) which has two kinds of expressions since two different separation schemes of 
running coupling and subtraction have been used in 0](see [is'] for more discussions on separation schemes). 
Balitsky took the transverse coordinate of either the quark at z±i or the antiquark at z±2 to be the subtraction point. 
He got the kernel of the running coupling contribution as follows @ 



i^«-Hr,ri,r2) 



27r2 



2 Jl 



rt r. 



1 '2 



1 f asjrl) 



1 



1 f asjr^) 



1 



(14) 



Here we introduce the notation r = x±_ — y±, ri = x± — z± and r2 = z± — yj_ for the sizes of parent and of the new 
daughter dipoles produced by the evolution. On the other hand, in the subtraction scheme proposed by Kovchegov- 
Weigert the subtraction point is fixed at the transverse coordinate of the gluon at z± — rjz^i + (1 — 'n)z±_2 in which rj is 
the longitudinal momentum fraction of gluon carried by quark. They got the modified kernel of the running coupling 
contribution 



K^^{V,V^,V2) 



27r2 



asirl)as{rj) ri ■ r2 



with 



i?^(r,ri,r2) = n r2 



^2 2 



(15) 



(16) 



The second term in r.h.s of (llip . S, which is referred to as 'subtraction' contribution, is given by 



S[S] 



(f z±i(f z±2 KQ{x±,y±; z±i, z±2) [S{xj_-wj_,Y) S{w±-y±,Y) - S{x±- z±i,Y) S{z±2-y±,Y)] (17) 



with the bare coupling. The interaction structures are modified in the above equation since the quark-antiquark 
pair is added to the evolved wave function (see Fig.2B). The JC(i){x±m, x±n; z±i, z±2) is a resummed JIMWLK kernel 
which can be found in 1131 



K(i){xj_,y±;z±i,z±2) = C'f ^ (-1)""+" /C(r)(a;±„i, a:;^„; 2^2)- 



(18) 



ni.n—O 
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In terms of Balitsky' subtraction scheme one substitutes w± = z^i or w± = zj_2 in Eq. ([T7)) and gets the subtraction 
term 

S^"-'-[S] ^ J (fz±i tfz±2 K(i){x±,y±; z±i,z±2) [S{x± - z±i, Y) S{z±i - y±,Y) - S{x± - z±i, Y) S{z±2 - y±, Y)]. 

(19) 

According to Kovchegov-Weigert's subtraction scheme one substitutes w± = z± — r]z±i + (1 — rj)^±2 in Eq. pT)) and 
gets 

S^'^[S] = J (fz±i (fz±2 K(i){x±,y±; z±i, z±2) [S{x± - z±,Y) S{z± - y±,Y) - S{x± - z±i,Y) S{z±2 - y±,Y)]. 

(20) 

Eq. (fT7|) shows that the S[S] is of order while Tl[S] is of order as and all terms of S[S] are quadratic in S, 
S{x± — w±,Y) S{w± — y±,Y), S{x± — z±i, Y) S{z±2 — y±,Y). Thus, for high rapidity and small S, the subtraction 
term is small as compared to the running coupling term, as also shown numerically in ^13 ] . Since this is the kinematic 
region in which we are interested in this paper, we hereafter will neglect the subtraction term. . In this paper we study 
the evolution equation in the saturation regime where the evolution equation including running coupling corrections 
can be solved analytically. 



B. Solution to Balitsky and Kovchegov-Weigert equations in the saturation regime 



In the saturation regime in which the interaction between partons is very strong, S{xj_ — y±, Y) — > 0, and unitarity 
corrections become important, the quadratic terms in (fTT|) can be neglected in which case one needs only keep the 
second term on the r.h.s of (jl2p . The evolution equation including running coupling is given by 



dS(x^^y^,Y) 
dY 



I 



(fz± K{r, ri, r2)5'(a;_L - y±,Y) 



(21) 



with modified kernel ii'(r, ri, r2). In the saturation region, rQs{Y) ^ 1, the main contribution to the integration on 
the r.h.s of (EU) comes from either 



1/Qs < ri < r; 



r2 r 



(22) 



1/Qs < r2 < r 



ri ^ r. 



(23) 



Let us look at one of them, i.e., when 1/Qs <C ri <C r, the r2 is approximate equal to r, r2 ~ r, the A'^'^'(r, ri, r2) 
becomes 



/^^"'(r,ri,r2) 



27r rf 



1 as{rl) 1 / as{r'^) 



rf as{r'^) ri \as{r() 



and the K^^(r, ri, r2) has the form as follows 



if^^(r,ri,r2) 



27r2 



. 2n 1 o "s(r^)a^(ri) ri ■ rz , 2n 1 



(24) 



(25) 



here we use i?^(r, ri, r2) ~ which can be obtained via simple calculation in (fT6|) with condition oi 1/Qs ri r 
and r2 ~ r. In the ri <C r2 limit it is the first term which dominates Eq. (|25p and has the running coupling scale 
given by the size of the smaller dipole 



i^^^(r,ri,r2) 



— ^a, ri — 

ZTT'' rf 



(26) 



We wish to note that the modified Balitsky and Kovchegov-Weigert kernels including running coupling have the same 
form in the saturation regime. It is an interesting outcome which means that the evolution equation with running 
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coupling corrections is independent of the choice of transverse coordinate of subtraction point in the saturation regime. 
And the modified BaHtsky and Kovchegov-Weigert equations with running couphng corrections are equivalent to each 
other in the saturation region. In other words, the S-matrix of the Balitsky and Kovchegov-Weigert equations are 
exactly the same in the saturation region. 

Now let us put the modified kernel (|24p or ((26l) into pT|) . we can get a simplified evolution equation as follows: 



dSjr, Y) _ Nc 
dY ~ 

with the running coupling at one loop accuracy 

«s(r?) = 



1 



d'n a,(rf)-5(r, Y), 



1 + Ml In ( ^ 



(27) 



(28) 



givmg 



dS{r, Y) _ N,^i 



dY 



TTfll 



In 1 + /^i In 



whose solution (see also [13]) is 



S{r, y) = e 
with the saturation momentum 



A2 



— In 1 + /ii In 



H — ^ ln( In ■ 



Ml 



S{r,Y) 



iS{ro,Y) 



In 



(29) 



(30) 



(31) 



We wish to note that the analytic result for the S'-matrix including the running coupling corrections is different as 
compared to the fixed coupling case. The exponent in Eq. (j30p is decreasing linearly with rapidity while the exponent 
in Eq. ([8|) is decreasing quadratically with rapidity, which indicates that the running coupling slows down the evolution 
of the scattering amplitude with rapidity. 



IV. EFFECTS OF RARE FLUCTUATIONS 



A. Fixed coupling case 

At very high energy the typical configuration of a dipolc's light-cone wavefunction is a Color Glass Condensate 
which is a state having high occupancy for all gluonic levels of momentum less than or equal to saturation momentum 
Qs- In the fixed coupling case, the authors of Ref. [5,] computed the S'-matrix of two typical configurations (of 
condensate type) and of dipole-typical configuration scattering, they found that the typical configurations lead to too 
small results for the 5-matrix, being proportional to exp{— ciQlrp/aj} and exp{— ^ {Q^Tq)} , respectively, ci and 
c are constant which are not important for our purpose. Thus they tried to search for configurations which are more 
rare in the wavefunction but which dominate very high energy dipole-dipole scattering and lead to a larger S'-matrix. 
They found the reason why the typical configurations have given a small S-matrix is that the typical configurations 
contain too many gluons at the time of collision, therefore leading to the S'-matrix is extremely small. This suggests 
that the strategy for finding the rare configuration is to minimize the number of gluons by suppressing the evolution 
(see next section for the details of how to obtain the rare configuration). The rare configuration is a state which have 
no more than one dipole of size nro or larger (with k a constant of order 1 and ro a size of parent dipole) when the 
system undertakes BK evolution. In the center of mass frame, the S-matrix is then given by the probability of the 
rare configurations for each of the parent dipoles partaking in the collision, times the S'-matrix for the scattering of 
two dipoles separated by a rapidity gap Yq, 

Sy «e-^i"'(0^o)5y„(ro) (32) 
which is significantly larger than the results coming from the condensate-condensate and dipole-condensate scattering. 
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Y/2- 



ik'rl)=\^c(Y-y /2) 







^ ln(k'rl)=--^c(y+yJ2) 



-Y/2- 



FIG. 3: The configuration in center of mass frame. 



B. Running coupling case 

Following the framework of Ref. [s*], consider the high-energy scattering of dipoles at zero impact parameter in 
the center of mass frame where one of dipoles is left-moving and the other is right-moving. In order to obtain rare 
configuration, we suppose that the wavefunction of the right-moving dipole consists only of the parent dipole with 
size ro in the rapidity interval Yq/2 < y < F/2, where Yq is the critical value of rapidity for the onset of unitarity 
corrections, with the similar requirement on the left-moving dipole in the rapidity interval —Yjl < y < —YqII. In 
the rapidity interval < y < Iq/S and — lo/2 < y < the right-moving and left-moving dipoles have normal BFKL 
evolution, respectively. 

However, we cann't require that all evolution of right-moving dipoles are absent in the rapidity interval lo/2 < 
y <Y . What we can do is to only allow that evolution which produces very small dipoles, in order to guarantee the 
system have no more than one dipole of size ktq or larger, with n a constant of order one. And we setup constraints to 
suppress the creation of dipoles much smaller than rg at rapidities y > Yq/S to avoid dipoles emitted at intermediate 
rapidities evolving into dipoles of size tq or larger at rapidity Y jl. We require that the gluon emission from the parent 
dipoles is forbidden if the gluon has k\_ and y in the shaded triangles of Fig. 3. The line 



ln(A;^ro2) = ^c(y-^) (33) 

and a similar line for the lower triangle, is determined by the requirement that giuons in the right hand side of that 
line cann't evolve by normal BFKL evolution into shaded triangles. 

Now we compute the probability of rare configurations §(j:_l — yj_,i^ ~ ^o) which has the same meaning as the 
survival probability of the parent dipoles after a BFKL evolution over a rapidities interval Y — Y^ Q . This probability 
decreases with increasing Y due to gluon emission and the corresponding rate is the same as the virtual term in (|13|) : 

^§(a;i -y^,Y-Yo) = - f d^z^ k{x^,v^, z^)^x^ ~ y±,Y - Yo) (34) 
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whose solution is similar to Eq. ([30] 



S{r,Y-Yo) 



In^ 



i+m 1, 



ln( 



(35) 



Let §(ro, (y — Yq)/2) denote the probability of a parent dipole not given rise to any emission of gluon in the upper 
triangle of Fig. 3. The S'-matrix can be obtained by the product of §(ro, {Y — Yq)/2) for each of the parent dipoles 
participating in the scattering, times S'(ro,lo) which is a S'-matrix for the scattering of two elementary dipoles. By 
using (pSjl. one gets: 



S{r,Y) = e 



In 



V2 



\/2 li 



■ In 



Siro,Ya) 



(36) 



which only brings in very small corrections to (j30p and indicates that the rare fluctuations are less important in 
the running coupling case as compared to the fixed coupling case 's'l, where the rare fluctuations are important and 
the exponential factor of S'-matrix in the saturation regime has twice as large as the result which emerges when 
fluctuations are taken into account. We also consider the rare fluctuations on top of the running coupling effects in a 
general frame (please see the Appendix), we find the same result as p6p . 



V. SHAPE OF DIPOLE CROSS SECTION WITH RUNNING COUPLING 

The authors of Ref. [l5| computed the scattering amplitude for rQs ^ 1 using BFKL evolution and running 
coupling. Combining the outcome of Ref. and our result (l30|) which is valid for rQs ^ 1, the shape of dipole 
cross section with running coupling reads: 



N{r,Y) 



M 



1 - e 



-(l-Ao) 



1-(S)+T4k 



Qs(y) 



rQs < 1 



So rQs > 1 , 



where Aq is the solution to x ('^o)(l — ''^o) = — x(-^o) with x the usual BFKL eigenvalue function, A^o is a constant but 
with no control at this moment and Qs is the saturation momentum including running coupling corrections. 
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APPENDIX A: RARE FLUCTUATIONS IN A GENERAL FRAME 



Consider a high energy scattering of a right-moving dipole of size tq and rapidity 1^ — 12 on a left-moving dipole 
of size ri and rapidity —Y2 in an arbitrary frame. The frame and scattering picture are illustrated in Fig. 4, where 
Yq is a rapidity gap between two dipoles. For later convenience, we require that Y2 < ■^(Y — Yq). We require that 
no additional dipoles can be created from the gluon emission of left-moving dipole ri which would have a strong 
interaction with the right-moving dipoles. The dipoles which would have such strong interactions would be of size 
r > l/Qs at the scattering time. So we should suppress the emission of those dipoles which could become of size 
1/Os or larger after a normal evolution over the rapidity interval —Y2 < y < 0. For the right-moving dipole ro, we 
suppress evolution over its Y — Y2 — (Yi + Yq) with the region of suppression given by the upper shaded triangle of 
Fig.4. The line 



Hk^rl) = ^c{y - Fi - Yq) 



(Al) 
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and a similar line for the lower triangle, is determined by the requirement that gluons locating in the right hand side 
of that line cann't evolve by normal BFKL evolution into shaded triangles. We will determine Yi by maximizing the 
5-matrix later. The unshaded triangle, whose rapidity values go from to Yl, is a saturation region where the dipole 
tq has evolved into a Color Glass Condensate. 

After we have a clear scattering picture of dipoles, the S'-matrix can be evaluated at hand 



5(ro, ri, F) = §fl(ro, F - Fo - ^ - Y2)S{ro,n,Yo + ri)Si(ri, ^2) 



(A2) 



with iS(ro,ri,lo + Yi) is the S'-matrix for scattering of a elementary dipole ri on a Color Glass Condensate state 
which is evolved from dipole vq and E>ii{ro,Y — Yq — Yi — Y2) and SL(ri,l2) are the suppression factor from the no 
emission requirement for two dipoles, which are given in terms of the area of the upper and lower shaded regions of 
Fig. 4. After using (??), one obtains 



SR{ro,Y-Yo-Y,-Y2) ^ e 
and 



(A3) 



§L(ri,y2) 



exp 



— cln 



cln — , . 1 . (ri + ^2) ^ 5- In 1 + ^1 Vc(yi + Y2) 



1 + 



Ml MI ^ 



1 + Mi\/cyi 



(A4) 



The S can be computed by using the BK equation with running coupling corrections since the BK equation with 
running coupling corrections describes correctly the scattering of an elementary dipole on a Color Glass Condensate. 
By using one gets 



i+Mi 1 



5(ro,ri,yo+n) = e 
Substituting (|A3]l . (|A4l) and ((A5|l into ((A2|l . one obtains: 

■l + MiVc(i^-^2-i"i->'o) 



5(ro,ri,y) 



exp 


iVcM 


cln 1^ 




CTTMl 





1 + Miln(;:4^) 



- 2 - - Fi - yo; 



S(ro,yo). (A5) 

, ^c{Y -Y2-Y^- Fo) 



Ml 



A In (l + MiVc(r-r2-Fi-ro)) +cln ^ + _ ^\ ^y^^y^^ + 

^in(i + Aiiv/c(ri + y2)) 



v/c(yi + F2) 



Ml 



S(ro,ro). 



(A6) 



which connects to a set of configurations of the wavefunction described by rapidity Y\ . The S'-matrix is determined by 
the values of Y\ which maximizes the r.h.s of Eq. \M>\ or equivalently minimizes the exponent of the r.h.s of Eq. (|A6|1 . 
We obtain 



Y^ = \{Y -Y,,)-Y2. 



(A7) 



Take this Y\ into \M>\ . finally the S-matrix is: 



In-" 



-k + 



S(r, r) = e 

which is exactly the same as the corresponding result 



^>n(l.^l„(4-)) 



(A8) 



in the center of mass frame. 
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